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Abstract. We reconsider the rejuvenation mechanism as proposed by 
Shull, Fesen, &: Saken (1989). These authors suggest that an active pulsar 
can catch up with, and rejuvenate the shell of the associated supernova 
remnant. The morphology of the SNRs G5.4-1.2 and GTB80 seem to 
confirm this rejuvenation mechanism. The spindown energy is deposited 
by the pulsar as a relativistic pulsar wind, and has a sufficient power to 
explain the observed radio emission observed in these remnants. Shull et 
al. (1989) did not explain the observed lengthscales of the rejuvenated 
parts of the SNR shell, therefore one needs to consider the diffusive 
transport of the injected electrons by the pulsar wind. We propose to ap¬ 
ply a diffusion mechanism as introduced by Jokipii (1987), which makes 
a distinction between diffusion along the magnetic field lines and per¬ 
pendicular to the magnetic field lines, parameterised by the gyro factor 
rj. We show that one has to assume a high value for the gyro factor, 
r/ ~ 10^ — 10^, i.e. diffusion of the electrons along the magnetic field line 
is much faster then perpendicular to the magnetic field line, in order for 
the rejuvenation mechanism to work on the observed lengthscales. 


1. Introduction 

A supernova remnant (SNR) results from the supernova explosion of a massive 
star. When the core of the progenitor star collapses, a neutron star can be 
formed. Several mechanisms proposed for a core collapse supernova can impart 
a kick velocity to the pulsar, although no particular mechanism can be favoured. 

The relativistic pulsar wind will interact with the SNR, resulting in a pulsar 
wind nebula (PWN), which contains the shocked pulsar wind material. Initially 
the PWN will be centrally located in the SNR, but due to the kick velocity of 
the pulsar, the PWN will be dragged along by the pulsar, being deformed into a 
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bow shock and ultimately break through the shell of the decelerating SNR (van 
der Swaluw, Achterberg, & Gallant 1998; Chevalier 1998). 

Figure 1 depicts the configuration of this last interacting stage of such a 
composite remnant. The relativistic pulsar wind is terminated by a strong MHD 
termination shock, whereas the bubble around the termination shock has been 
deformed into a bow shock due to the supersonic motion of the pulsar. Rel¬ 
ativistic particles are injected at the site of the termination shock, of which a 
large fraction is advected away resulting in a wake of relativistic particles, which 
explains the observed trail of radio emission for bow shock PWN. However, due 
to the diffusive transport of the injected particles, a small fraction of these parti¬ 
cles will radiate part of their energy away in the magnetic field of the SNR shell. 
Shull et al. (1989) have argued that these freshly injected particles brighten the 
radio emission of the SNR shell. 

We calculate the diffusive length scales of the injected electrons in the life¬ 
time of a SNR or their radiative lifetime. This will enable us to compare these 
results with actual observed lengthscales in these systems. We will show that 
we need strongly anisotropic diffusion in order to obtain agreement between our 
calculations and the observed properties of several composite remnants. 

2. Synchrotron theory and diffusive transport 
2.1. General theory 

Relativistic electrons in an astrophysical flow radiate part of their energy away as 
synchrotron radiation, due to the presence of a magnetic field, B in the plasma. 
Using standard synchrotron theory (see e.g. Rybicky & Lightman 1979), one 
can write the frequency where the emission peaks as: 

z^MHz — 4.67R^G-EGeV MHz, (1) 

here F'ceV is the energy of the electron in GeV. The timescale on which the 
electron has lost half of its energy due to synchrotron losses can be written as: 

Tioss - 2 X yr. (2) 

The propagation of relativistic particles through the flow of a plasma is a com¬ 
bination of advection by the large-scale flow, and diffusion with respect to this 
flow. We will consider the limit of Bohm diffusion, in which case the mean free 
path A equals the gyroradius rg of the particle. Using this limit, one can rewrite 
the Bohm diffusion coefficient kb in terms of the above characteristic frequency 
z^MHz, which yields: 

Kb - 7.5 X cm^/sec. (3) 

Now by considering particles injected at t = 0, one can write the diffusion length- 
scale Aiigyn, at the synchrotron loss time rioss in the case for Bohm diffusion 
as: 

ARsyn = \/2 Kb rioss - 9.7 R“q^^ parsec. (4) 

The above equation will serve as a reference value, when we consider the more 
general case of diffusion in the next subsection, where the mean free path A 
satisfies A » rg. 
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Figure 1. Configuration of a pulsar wind interacting with the shell 
of a supernova remnant. The pulsar wind is terminated by a strong 
MHD shock (dashed line) and the PWN itself is bounded by a bow 
shock (solid line). The PWN is propagating through the shell of its 
associated SNR, where the magnetic field lines are parallel with the 
shell of the SNR. By considering anisotropic diffusion, one cane explain 
the lengthscales of rejuvenated shells of SNRs like G5.4-1.2 and CTB80. 


2.2. Application to pulsar winds 

We consider the case of a pulsar wind close to the shell of a SNR. Because of the 
small size of a PWN bow shock (~ 0.1 parsec) compared with the SNR (~ 10.0 
parsec), we can approximate the site where the relativistic electrons are injected 
as a point source in the SNR interior. We can calculate the diffusion lengthscale 
again, this time using the crossing time for a pulsar with speed Rpsr to catch up 
with the SNR shell with radius Rsm, he. Rsnr/Fpsr. Using Bohm diffusion again 
yields; 

ARradio - 2.2 X 10“^jz^^H^R^Q/^(Rsnr/Upsr)^/^ parsec. (5) 

From the above equation one can see, that using standard Bohm diffusion, the 
obtained dijfusive lengthseale is similar to the size of the PWN itself i.e. the 
radio electrons are confined to the surroundings of the PWN. However, the 
observed size of the rejuvenated shells is much larger (~ 10 parsec), therefore 
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we have to extend the diffusion model in order for the diffusive lengthscales to 
match the observed lengthscales. 

We follow Jokipii (1987) by making a distinction between diffusion along 
the magnetic field lines (k||) and perpendicular to the magnetic field lines (k_l). 
The mean free path along the magnetic field line then equals Ay = rjrg, with rj 
the gyro factor, which is related to the turbulence level 6B in the magnetic field 
by Tj = {6B/B)~^. The diffusion coefficient along the magnetic field lines equals 
Ky = r]KB- One usually assumes 77 ^ 1. Perpendicular to the magnetic field lines 
we follow Jokipii (1987) and assume that a particle scatters one gyroradius across 
field lines for every parallel scattering length, so k_l = 77Kb/(1 + f]'^) — kb/ti- 
With this description for diffusion, particles diffuse faster along the magnetic 
field lines then compared with quasi-isotropic Bohm diffusion. Using this de¬ 
scription together with the timescale for synchrotron losses, Tiogg, one can write 
the diffusion lengthscale at time t of the radio electrons as: 

Ai?y = AR,yn X ^ X - , ( 6 ) 

\kbJ \nossJ 

with rioss given by equation (^) and AE^yn given by equation (Q). The above 
equation gives the diffusive lengthscale along the magnetic field, where the dif¬ 
fusion proceeds rapidly. 

2.3. Comparison with observations 

In this section we will use equation (^) to determine lower limits for the gyro 
factor 7] by matching Ai?|| with the size of the rejuvenated shells observed in 
the SNRs CTB80 and G5.4-1.2. We make two assumptions which will lead to 
a minimum value for the gyro factor: 1 ) the timescale t for the interaction 
between the pulsar wind and the shell of the remnant is taken to be the age 
of the remnant; 2 ) the diffusion process is assumed to take place in a uniform 
magnetic field rather then in a curved, position dependent magnetic field. 

CTB80 and PSR 1951+32 We use as a reference Strom k. Stappers (2000), 
who observed this remnant at a frequency of = 600 MHz. Using a distance 
of 2 kpc, the rejuvenated shell has a size of ~ 17 parsec. Taking the age of the 
system as t = 10^ years, we obtain rioss — 8.2 x B^q"^ years. Using equation 
(|^, we derive a lower limit for the gyro factor in order to match this observation: 

77 > 2.5 X liy^B^J^. (7) 

G5.4-1-2 and PSR B1757-24 We use Gaensler & Frail (2000) as a reference 
for this remnant. Their map is at a frequency of 330 MHz. The rejuvenated 
shell of the SNR has a size of ~ 30 parsec for a distance of 4.0 kpc. By using a 
characteristic age of t = 1.6 x 10 “^ year, we derive a value rioss — 1-1 x 
years. Using equation (I) again, we derive a lower limit for the gyro factor, 
which equals: 


77 > 6.6 X lO^Rjg- 


(8) 
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Figure 2. Synchrotron maps illustrating the effect of different angles 
of observation with a gyro factor rj = 10. 


3. Simulations 

In this section we show results from a Monte Carlo simulation, which traces 
the propagation of test particles in the flow of a bow shock. We obtained the 
steady-state flow inside a bow shock by performing hydrodynamical simulations 
with the Versatile Advection Code [](T6th 1996). 

The test particles are contineously injected in the flow of this bow shock at 
the site of the termination shock. We use Ito stochastic differential equations 
(SDEs), to simulate the random walk trajectories of these particles in the flow 
of the bow shock (a detailed discussion can be found in the Ph.D. thesis of E. 
van der Swaluw 

By considering many realizations of the SDE in the bow shock flow, we 
obtain a steady-state distribution of particles in phase space. We use this dis- 


^See ^ittp:/ /www.phys.uu.nl/~toth/ 


publicly available at: http://www.library.uu.nl/digiarchief/dip/diss/1967493/inhoud.litm 
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tribution to produce a synchrotron map, which can illustrate the influence of a 
gyro factor rj > 1. 

The hydrodynamical flow is axially symmetric around the Z-axis and the 
magnetic field strength is uniform inside the PWN. Outside the PWN, the mag¬ 
netic field is uniform and directed perpendicular to the pulsar velocity (see also 
figure 1). 

Figure 2 shows the result for a case with r] = 10. In the upper left panel the 
line of sight of the observer is perpendicular with the magnetic field lines outside 
the PWN. In the lower right panel the line of sight is parallel with the magnetic 
field lines, therefore the rejuvenated parts of the remnant are not visible. 

4. Conclusions 

We have reconsidered the rejuvenation mechanism as proposed by Shull et al. 
(1989). This was done by investigating the propagation of the injected rela¬ 
tivistic electrons at the site of the wind termination shock through the PWN 
and the associated SNR shell. We conclude that a pulsar wind can rejuvenate 
the shell of a SNR if the diffusive transport of the injected electrons is strongly 
anisotropic, i.e. the gyro factor has to have a minimum value of rj lO^ - 10^. 
The limits for the gyro factor derived in this work are lower limits because of the 
assumptions we made about the interaction time and the configuration of the 
magnetic field (section 2). Furthermore, our numerical experiment has shown 
that the line of sight between the observer and the magnetic field introduce a 
projection effect, which diminishes the lengthscales of the rejuvenated parts of 
the SNR shell (section 3). 
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